Quark distributions in the nucleon based on a 
relativistic 3-body approach to the NJL model 

H. Mineo, W. Bentz, K. Yazaki 
Department of Physics, 
Faculty of Science 
University of Tokyo 
Hongo 7-3-1, Bunkyo-ku, Tokyo 113, Japan 

Abstract 

Quark light cone momentum distributions in the nucleon are cal- 
culated in a relativistic 3-body approach to the NJL model by using a 
simple 'static approximation' for the Faddeev kernel. A method is pre- 
sented which automatically satisfies the number and momentum sum 
rules, even in the regularized theory. In order to assess the sensitivity 
to the regularization scheme, two schemes which can be formulated in 
terms of light cone variables are discussed. The effects of the (com- 
posite) pion cloud are taken into account in a convolution approach, 
and the violation of the Gottfried sum rule is discussed. After per- 
forming the Q 2 evolution, the resulting distributions are compared to 
the empirical ones. 
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1 Introduction 



Deep inelastic lepton-nucleon scattering experiments are providing detailed 
information on the quark and gluon distributions in the nucleon By 
analyzing the distribution functions derived from the measured cross sections 
with the help of the Q 2 evolution based on perturbative QCD, one can extract 
the fraction of the nucleon's momentum and spin carried by the quarks at 
some renormalization scale /i, and obtain valuable information on the spin- 
flavor structure of the nucleon [0]. For example, of great interest in this 
connection are the flavor dependences of the valence quark distributions 
and of the sea quark distributions j|, [5] , which reflect the roles of qq (diquark) 
and qq (mesonic) correlations in the nucleon. Current experiments at HERA 
are exploring the region of small x in order to get more information on large- 
distance (non-perturbative) phenomena, and future experiments at RHIC 
will concentrate on the spin distributions. 

An essential tool to analyze the data on the nucleon structure functions 
in the Bjorken limit are the factorization theorems ||, which allow us to sep- 
arate the long-distance parts (parton distributions) from the short-distance 
parts (hard scattering cross sections). The latter ones are calculable from 
perturbative QCD, while the parton distributions require the knowledge of 
the nucleon wave function. The central point of the factorization theorems 
is that the infrared- divergent (long-distance) contributions to the perturba- 
tive diagrams can be absorbed into the definition of 'renormalized' parton 
distributions @. Due to this procedure, the parton distributions eventually 
depend on a factorization scale (/i), which is usually taken to be the same 
as the renormalization scale. The requirement that the structure functions 
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should not depend on this scale then leads to the famous DGLAP equations 
[0] for the /i dependence of the distribution functions, which in turn also 
determines the Q 2 dependence of the structure functions. One can therefore 
use effective quark theories to calculate the distribution functions at some 
'low energy scale' \x = Qo, where a description in terms of quark degrees of 
freedom alone is expected to be valid, and then use the DGLAP equations to 
relate them to the distribution functions extracted from experimental data 
at fi = Q. 

As an effective quark theory in the low energy region, the Nambu - Jona- 
Lasinio (NJL) model is a powerful tool to investigate the properties of 
hadrons ||. It exhibits the spontaneous breaking of chiral symmetry in a 
simple and clear way, and allows the solution of the relativistic two- and 
three-body equations in the ladder approximation due to the simplicity of 
the interaction |10|, . Concerning the ground state properties of baryons, 
the Faddeev approach to the NJL model has been quite successful jT2|, and 
since in this approach the translational invariance and covariance are pre- 
served, it is natural to apply it also to the structure functions of the nu- 
cleon.[| Actually, the NJL model has been formulated recently on the light 



cone (LC) and the structure function of the pion has been calculated [[15 



For the case of the nucleon, however, there is the problem that the relativistic 
Faddeev equations are usually solved by performing a Wick rotation using 
the Euclidean sharp cut-off [JIDI , and this method cannot be applied directly 
to the calculation of the LC momentum distributions. A possibility is first 
to go to the moment space and then re-construct the distribution functions 



2 There are already several investigations of the nucleon structure functions using the 
soliton approach to the NJL model [fL3|. For a calculation using the Bethe-Salpeter equa- 
tion for a quark and a structureless diquark, see ref. mM. 



3 



by the inverse Mellin transformation. While such a calculation is now in 



progress [[16] , it is desirable for a first orientation to use some simple approx- 
imation which allows a direct calculation of the distributions. The 'static 
approximation' to the Faddeev kernel [T7], |T8|] , which amounts to neglect the 
momentum dependence of the propagator of the exchanged quark, allows an 



analytic solution of the Faddeev equation, and it has been shown in ref. [18 



that for the nucleon mass this approximation is not unreasonable. However, 
the overbinding of the nucleon compared to the exact Faddeev result and the 
'point-like' quark-diquark interaction lead to radii which are too small [|H| . 
and therefore we can expect that the resulting momentum distributions will 
be too stiff. However, our main purpose here is to see some general trends, 
to investigate the influence of the regularization scheme and to demonstrate 
that in the Faddeev approach the validity of the number and momentum sum 
rules is guaranteed from the outset, which is not the case in the soliton or 
bag model approaches. The calculation employing the static approximation 
to the Faddeev kernel resembles the quark-diquark calculations performed 
earlier 0]. We will go beyond these calculations by including the struc- 
ture of the diquark (and also of the pion when estimating the pionic cloud 
effects), and investigating the sensitivity to the regularization scheme while 
preserving the number and momentum sum rules. 

The rest of this paper is organized as follows: In sect. 2 we explain the 
model for the nucleon wave function, in sect. 3 we explain our method to 
calculate the quark distribution functions, and in sect. 4 we discuss the 
numerical results. A summary is presented in sect. 5. 
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2 NJL model for the nucleon wave function 



The NJL model is characterized by a chirally symmetric four- Fermi interac- 
tion lagrangian £j. By means of Fierz transformations, one can rewrite any 
Lj into a form where the interaction strength in a particular qq or qq channel 
can be read off directly [M . That part which generates the constituent quark 



mass Mq and the pion as a collective qq bound state is given by 

= ^((^) 2 - 0%5T)V>) 2 ), (2.1) 

and that which describes the qq interaction in the scalar diquark (J n = 
+ , T = 0) channel is 

d, s = 9s (TsC) T 2 p A i?) (> T (C-^s) r 2 p^) , (2.2) 



where /3 = y 3/2 A (A = 2, 5, 7) are the color 3 matrices, and C = 27270- 
The coupling constants g n and g s are related to the ones appearing in the 
original Li by numerical factors due to the Fierz transformation, but instead 
of choosing a particular form of Ci we will treat g n and the ratio r s = g s / 'g v 
as free parameters, where the latter reflects different possible forms of Li 



m 



The reduced t-matrices in the pionic and scalar diquark channels are 
obtained from the respective Bethe-Salpeter (BS) equations as |L0[ 



1 + 25- 7r n 7r (fc) 1 + 2g s U s (k) 



with the bubble graph 



l5S(qhbS(q - k) 



(2.4) 
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where S(q) = — is the Feynman propagator and Mn the con- 

4- Mq + 16 

stituent quark mass. 

If the interacting two-body channels are restricted to the scalar diquark 
one, the relativistic Faddeev equation |L0| can be reduced to an effective 



BS equation for a composite scalar diquark and a quark, interacting via 
quark exchange [2~I|. As we explained in the previous section, in this paper 



we will restrict ourselves to the static approximation |17|, [Tj|, where the 
Feynman propagator in the quark exchange kernel is simply replaced by 
—1/Mq. Then the effective BS equation reduces to a geometric series of 
quark-diquark bubble graphs (Hn(p)), and the solution for the t-matrix in 
the color singlet channel is 

™ = iTT^) < 2 ' 5 > 

with 

nN<p) = -fj^S{k)T.(p-k). (2.6) 

The quark-diquark vertex function Tn{p) in the covariant normalization is 
then obtained from the pole behaviour T — > J2 S ^n(ps)Tn(ps) / (p 2 —M N +ie), 
where Mn is the nucleon mass, as 

T N (p) = Z N u N {p) (2.7) 

Zn = = {u n (p)4^Mp)) ' (2 ' 8) 

where un{p) is a free Dirac spinor with mass normalized by un{p) un(p) = 

3 Our conventions for LC variables are o ± = -j= (a ± a 3 ) , a± = A= (ao ± 0.3) , and 
a\ = —a±i (i = 1, 2). The Lorentz scalar product is a ■ b = a + b + + a—b~ — a± ■ b±. We 
will frequently call p-(p+) the 'LC minus (plus) component' of the four vector p. 
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2Mn. In this normalization, the vertex function satisfies the relation 

1 r„( P )« r „ w = lf (2 . 9) 



<9p + 

which leads to charge and baryon number conservation in any treatment 
which preserves the Ward identity 

K, IP (P-P) = ^ JV./P (2.10) 

for the vertex of the quark number current of the proton at q = {N u /p = 
2, N d/P = 1). 



3 Quark distribution functions 

The twist- 2 quark LC momentum distribution in the proton (momentum p) 



is defined as [22, 23 



} ' q/p{x) = \S ^r e ">' T (M Q h + M*-)) \p), (3-i) 

where q denotes the quark flavor, \p) is the proton state, and x is the Bjorken 
variable which corresponds to the fraction of the proton's LC momentum 
component p_ carried by the quark q. As has been discussed in detail in ref. 



[p2] , |23| , for connected LC correlation functions the T-product is identical to 
the usual product, from which it follows that the distribution ( |3.1| ) is non-zero 
in the interval — 1 < x < 1. The physical quark and antiquark distributions 
which determine the structure functions Fi and F 2 are obtained for < 
x < 1 as f q /p(x) = f q /p(x) and fq/p(x) = —f q /p(—x). The valence (v) and 
sea (s) quark distributions are then given by f qv /p(x) = f q /p(x) — f q /p(x), 

4 This vertex is defined by (p|^(0)7 + i±^V(0)|p) = T N {p)A+ /p {p,p)T N {p) for q = u{d). 



The evaluation of the distribution ( |3.1| ) can be reduced to a straight 
forward Feynman diagram calculation by noting that it can be expressed as 

[0,11,110 

hipl?) = ^ ^ - y) Tr ° 1 + M q (p, k), (3.2) 

with the quark 2-point function in the proton given by 

M qtfia (p, k)=ij d 4 z e ikz < p\T $ q , a {0)il> q M) \P > ■ (3-3) 

We therefore have to evaluate the Feynman diagrams for the quark propaga- 
tor in the nucleon, trace it with 7 + , fix the LC minus component of the quark 
momentum as k- = p-X and integrate over the remaining components k + 
and k±. Since in our model for the nucleon discussed in the previous section 
the quark can either appear as a spectator or as a constituent of the scalar 
diquark, the Feynman diagrams to be evaluated are shown in fig.l. (In the 
full Faddeev approach, there is also a diagram where the external operator 
acts on the exchanged quark, but in the present static approximation this 
diagram does not contribute.) 

To present the formulae for the diagrams fig.l, we note that the sec- 
ond diagram (the diquark contribution) can be expressed conveniently as a 
convolution integral if we insert the identity 

1 = JdyJdzJ dql 5{y - ^) S(z - ^) 5{q 2 - q 2 ), (3.4) 

i.e; y is the fraction of the nucleon's momentum component p_ carried by the 
diquark, z is the fraction of the diquark's momentum component g_ carried 

5 The original Lorentz invariant expression is recovered by — > k ■ q/p ■ q and 

7 + /p_ — > tf/p ■ q. The expression ( |3.2| ) corresponds to the Bjorken limit in the frame 
where q ± =0, q z < 0, i.e; q+ — > oo, g_ — > —p_x. 
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by the quark inside the diquark (x = yz), and q% is the virtuality of the 
diquark. Using also the identity 

S(k)^S(k) = (3.5) 

and performing partial integrations in the plus components of the loop mo- 
menta, which is permissible since these integrations are convergent and not 
restricted by the regularization schemes to be discussed later, we obtain the 
following expression 6 : 

/q/pO) = 8 Q ,v F q/p {x) + ^F Q{D yp(x). (3.6) 

Here the first term corresponds to the first diagram in fig.l and is expressed 
as 



f Q/p( x ) = —F N ( — U N {x, p) ) T N , (3.7) 



-r N — 

2p_ \ dp + 

where Un(x,p) is the quark-diquark bubble graph for fixed minus component 
of the quark momentum: 

U N (x,p) = - J -^-5(x - y)S{k)r s {p - k). (3.8) 

The second term in (|3.6|) corresponding to the second diagram in fig.l is 
purely isoscalar and given by the convolution integral 

Fq(d)/p(x)= dy dzd(x-yz) dq 2 F Q/D (z,ql) F D / P (y,<fi), (3.9) 

JO JO J—oo 

where the distributions Fq/d(z, <2o) an d Fr>/p(y, g 2 ) for fixed virtuality of the 
diquark (D) are expressed as 

F Q/D (z,ql) = -2, 2 (g 2 )^il (3 . 10 ) 



6 To distinguish the case without the pion cloud (valence quark picture) from that 
including the pion cloud, we replace q — > Q in the formulae corresponding to the diagrams 
of fig.l. 



F D/p{y,ql) = T N (^—^- + y^ 1 )U N (y,q 2 ;p) T N . (3.11) 



2p_ dp + dq { 



Here H s (z, q$) is the quark-quark bubble graph for fixed minus momentum 
component of the quark, g 2 (q$) = — 1/ — ^ 2 ° J is the quark- diquark cou- 
pling constant, and Hjsr(y,qQ;p) is the quark-diquark bubble graph for fixed 



virtuality and minus momentum component of the diquark: 



^s{z,ql) 



d 4 k k_ 
(2tt) 4 g_ 

d 4 q q 



(3.12) 



-J q 2 =ql 



^N(y,q 2 ; P ) = -/^%-^)% 2 -g 2 )^-^(g). (3.13) 

If we use the dispersion representation for the diquark t-matrix r s , we can 
perform the k + and q + integrals analytically and verify that the distribution 
( |3.6|) is non-zero only in the interval < x < 1. This fact was anticipated 
already in our notation (/q/p(x) = /q/p(x)) and in the integration limits in 
eq. ( p.9|) , and corresponds to a valence quark model {fq/p{x) = 0). 

Using the above expressions and the normalization ( |2.9|) , it is a simple 
matter to confirm the validity of the number and momentum sum rules f\ 

f dx f Q/P (x) = N Q/P (3.14) 

f dx-x (fu/p(x) + f D /p(x)) = 1. (3.15) 

In more general terms, what we have really confirmed here is the validity 
of the Ward identities for quark number and momentum conservation: Con- 
cerning the number conservation, if we integrate the distribution function of 

7 To verify these relations, it is only necessary to note that the relations 

/ dzF Q/D (z,ql) = 2 and / dz ■ zF Q / D (z,<$) = 1 hold for any ql (here Fq/ d (z,$) is 

Jo Jo 

symmetric around z = 1/2), and therefore the second term in eq. (3.11) gives a vanishing 

surface term when integrated over The sum rules are then obvious since the integral 

of njv(y, <?o,2?) over q$ reduces to I1jv(1 — y,p), and that of IIjv(x, p) over x to Hn(p). 
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fig. 1 over x, the restriction = p_x is lifted, and the diagrams correspond 
to AQ/ p (p,p)/2p_. The validity of the Ward identity ( 2.10Q follows then from 



(3.5) and partial integrations in the plus components of the loop momenta. A 
similar argument holds for the Ward identity expressing momentum conser- 
vation. Therefore, the Ward identities and the sum rules (|3.14 ), ( 3.15Q hold 
in any regularization scheme which does not restrict the LC plus components 
of the loop momenta. The regularization schemes to be discussed at the end 
of this section satisfy this requirement. 

As we have noted above, the model described so far gives essentially 
only valence-like distributions at the low energy scale. Although sea quark 
distributions will be generated in the process of the Q 2 evolution, those will 
be flavor independent (fu/p = fd/p)i which contradicts the experimentally 
measured violation of the Gottfried sum rule JIJ. Also, it is a general trend 
of valence quark models that the resulting valence quark distributions are 
too stiff (too strongly pronounced peak and too small support at low values 
of x). We therefore consider here the effects of the pion dressing of the 



constituent quarks, as has been done also in previous works [25]. (For a 
recent investigation using the LC quantization, see ref. [2H|.) In order to 



take into account the pion cloud, in principle we should solve the Schwinger- 
Dyson equation for the quark Feynman propagator S(q) = — — , 

4- Mq-e q {p) 

where 

S Q (p) = -3 J-0- (75^)75) UV -k) = -Uq(p) (3.16) 

is the quark self energy due to the pion cloud, and the reduced pion t-matrix 
f n = r 7T + 2ig 7r depends also on S(p). This propagator should then be used to 
calculate the diquark t-matrix r s and the nucleon vertex function T^- Using 
these modified propagators and vertex functions, one should, in addition to 
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the diagrams of fig.l, also evaluate the diagrams of fig. 2, where the operator 
insertion is made on the quark while the pion is 'in flight', or on the quark 
and the antiquark in the pion. 

Clearly, such a calculation is very complicated, and the usually employed 



convolution formalism |25[ involves the following two major approximations: 
First, pionic effects can be renormalized into a re-definition of the constituent 
quark mass and the four fermi coupling constants if one approximates the 
quark propagator by its pole part Q 

S(p) = l Q , = Z Q S(p) (3.17) 



with 

z «=( 1+r i^j • < 3 - 18) 

If we define 'renormalized' coupling constants G a by G a = G /Zq (a = 

n, s), it is easy to see that the Green functions and vertex functions are 

renormalized according to r Q = T a /Zg, T = Zq T,T n = ^J~ZqTn- If we 

then impose the same conditions on the parameters as commonly used in 

the case without pion cloud effects (that is, /„- = 93MeV, m n = 140MeV, 

M Q = 300 - 500Me1/, and M N = 9A0MeV), the cut-off and the coupling 

constants G a take the same values as in the case without pionic cloud effects. 

(In terms of the lagrangian, such a renormalization procedure corresponds 

to writing G (V?/>) — G (ipipj with ip = ^J~Zq ijj and G = G/Zq.) It is also 

easy to check that all diagrams in fig.l and fig. 2 get a factor Zq. 

8 Green functions and vertex functions which differ from those without pionic effects 
by the replacements Mq — > Mq, G a — > G a = G a Zg (a — w,s) will be denoted by 
a hat. Since due to the discussion following eq.( 3.18 ) these renormalized quantities are 
numerically equivalent to the ones used previously, this distinction will eventually be 
dropped. 
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If one then calculates the diagrams of fig. 2 in terms of these renormalized 
quantities and writes the results in terms of a convolution integral, one finds 
that due to the Dirac structure of the insertions on the 'parent' quark line 



off shell case (fQ/p(x) — > fQ/p(x, k^), where k$ is the virtuality of the parent 
quark), convoluted with the quark distribution within the parent quark. Each 
of the other two terms involve one additional factor of (k$ — Mq) in the 
integrand compared to the first term, which has a sharp peak at k% = Mq. 
The second approximation commonly used is therefore to neglect these two 
terms, and to assume that due to the sharp peak of /q/p(x, /cq) the quark 
distribution within the parent quark can be evaluated at k\ = Mq and 
taken outside of the k% integral [p2]| . In this way one arrives at the familiar 
convolution form 



and a similar expression with q — > q, where the parent quark distribution in 
the proton f Q / P is given by eq. (|5T6|) , and f q /Q {fq/o) is the quark (antiquark) 
distribution within an on-shell parent quark which is obtained by evaluating 
the Feynman diagrams shown in fig. 3. 

The quark and antiquark distributions in the parent quark obtained from 
the diagrams of fig. 3 can be expressed as 



there appear three convolutions terms [p7[] . Only one of them involves the 
generalization of the 'bare' quark distribution in the nucleon (eq.( |3.6|) ) to the 




Q=U,D 



(3.19) 




(3.20) 



fu/u(x) 



fd/u(x) 



-F q/Q (x) + -F q{w)/Q (x) 
1 



(3.21) 



(3.22) 
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ft/u( x ) = Q F g(n)/Q(x). (3.23) 

The distributions for Q = D are also determined from these expressions due 
to isospin symmetry. The detailed formulae for the distributions F q /q(x) and 

rl rl roo 

^M/qO) = / dy dz5(x-yz) dq 2 Q F q /^{z,ql) F„ /Q (y,ql), (3.24) 

JO JO J—oo 

corresponding to the second and third diagrams of fig. 3, respectively, can be 
obtained from the previous expressions Q3.7D -( [3T3l) as follows: F q /Q(x), F q / n 
and F w /q are given by expressions similar to eqs. (|3.7|) , ( [3.10| ) and (|3.11|) , 
respectively, but with the following replacements: (i) The nucleon spinor T N 
is replaced by the quark spinor defined via the residue of the propagator 
S(p): Tq(p) = Zq Uq(p), where uq(p) is a free Dirac spinor with mass Mq 
normalized by uq{p) uq(p) = 2Mq. From eq.( [3.18| ) it follows that this spinor 
satisfies the relation 

1 r Q ( P )^M r o(p) = i - z Q . (3.25) 



2p_ dp+ 

(ii) the polarizations Ujsr(x,p) and IIiv(y,?ojP) are replaced by Uq(x,p) and 
Hq(?/, q^iP). These are defined analogously to eqs. fl3-8|) and ( p,13|) by intro- 
ducing the 5 function insertions to fix the minus momentum components of 
the quark or the pion and the virtuality q^ of the pion into Hq defined by 
eq. (|3.16l ) instead of II jv, and (iii) F g / n is given by the r.h.s. of eq. (|3.10|) , 
but without the overall factor 2. We note that f q /Q, and therefore also the 
distribution f q / P of eq. fl3.19|) , involves an overall factor Zq, in accordance 
with our discussion following eq. Q3.18Q . 

^From these expressions and the normalization (3.25), the validity of the 
number and momentum sum rules 

/ dx (/,/o(ar) - fq/g{x)) = 6 Q>q (3.26) 

J 
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dx ■ x (f u / Q {x) + fu/ Q {x) + fd/ Q {x) + fd/ Q {x)) = 1 (3.27) 



can be easily checked in the same way as outlined above for the parent 
quark distributions in the proton. (The number sum rule is a consequence 
of the Ward identity for the quark number current of the parent quark 
A^q(p,p) = — ^ ^ 5 qj Q.) The validity of the number and momentum sum 
rules 

/ dx f gv /p(x) = f dx (f q /p{x) - fq/p(x)) = N q/P (3.28) 

J «■> 

f dx-x (f u /p{x) + h/p{x) + f d /p{x) + fz/p( x )) = 1 ( 3 -29) 

J 



is then a consequence of eqs. ( |3.14| ), ( |3.15| ), ( |3.26| ) , refbmom2 and ( p,19| ). 
Of particular interest is also the Gottfried sum 

Sg = \j Q dx (f u /p(x) + fu /P {x) - fd/p{x) - fd/ P (x] 

= l --lf o dxF q/Q {x) = l --^{l-Z Q ), (3.30) 

which shows that the deviation from the valence quark model result (Sq = |) 
is due to the decrease of the probability of the 'bare' valence quark state 



:z Q <i) m. 



We now discuss our regularization scheme. Since the above expressions for 
the quark distributions involve loop integrals with one of the LC momentum 
components fixed, it is clear that we need a regularization scheme which can 
be formulated in terms of LC momenta. Two such schemes which have been 



discussed extensively in ref. [15| are the Lepage-Brodsky (LB) or invariant 
mass scheme , and the transverse cut-off (TR) scheme 137J. The basic 
graphs which are regularized in both schemes are the qq and qq bubble graphs 
Il s = ![„-, the quark-diquark bubble graph LT^r and the quark self energy LTq, 

15 



either for the case that all internal momentum components are integrated 
out or one of the LC momentum components is fixed. Concerning the LB 
scheme, it has been shown in detail for the case of U s in ref. [[15[ that, if 
all momentum components are integrated out, this scheme is equivalent to 
the covariant 3-momentum (or dispersion) cut-off scheme. f\ Generally, if the 
intermediate state involves particles with masses mi and rri2, the LB cut-off 
applied in the frame where the transverse components of the total momentum 
are zero (pj_ = 0) restricts the invariant mass of the state according to 

ki±^ + kI±£i <AiB> (3 . 31) 

where x and 1 — x are the fractions of the total momentum component p- 
carried by the two particles. The LB regulator Alb is related to the 3- 



momentum cut-off A 3 by Alb = y\/ m i + + ym^ + A3 J. For the case of 

Il s = U n we have mi = m 2 = Mq, and the value of A3 is determined as usual 

by requiring that /„- = 93MeV. In the case of the graphs ILv (IIq) we have 

rrii — Mq, while m 2 is the mass parameter in the dispersion representation 

of t s (tv). In order not to increase the number of parameters, we will take 

the same value of A3 for all graphs II S = IL, ILv and ILq. 

Concerning the TR cut-off scheme, it has been discussed in ref. JOJ that 

the use of this scheme requires a mass renormalization procedure, since the 

basic self energies IT a (a = s, tt, N, Q) involve also logarithmically divergent 

longitudinal momentum (k_) integrals, which are not affected by the TR 

9 By 'covariant 3-momentum cut-off scheme' we mean the procedure where the 3- 
momentum cut-off is introduced in the particular Lorentz frame where the total momentum 
of the two-body (qq, qq, or quark-diquark) state is zero (p = 0), and the result is 'boosted' 
to a general frame. For the graphs II Q (a — s, tt, N, Q) considered here, this 'boosting' 
simply means the replacement Pq — > p 2 , and for I±at and IIq also po7° — * i>- It is known 
that this procedure is equivalent to the dispersion cut-off scheme Q. 
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regularization prescription |fej_| < Atr- In this scheme one has therefore 
to impose the pion, the diquark and the nucleon masses as renormalization 
conditions rather than to relate them to the parameters g n , g s and r s via the 
pole conditions. For example, if we impose the condition 1 + 2g s II s (M D ) = 



for some fixed Md, the t-matrix r s in eq. (|2.3| ) can be re- written in the 
renormalized form r s = 2i / (U s (k 2 ) — n s (M|,)), which is formally independent 
of g s and free of divergencies due to the longitudinal momentum integration. 
In the calculation using the TR cut-off we will impose the same value of Md 
as obtained in the calculation using the LB cut-off. 

4 Numerical results 

In both the LB and TR regularization schemes we use Mq = 400MeV for the 
constituent quark mass, and determine the cut-off so as to reproduce f n = 
93MeV. This gives A 3 = 593MeV for the equivalent 3-momentum cut-off in 
the LB scheme, and Atr, = 407 MeV in the TR scheme. In the LB scheme, we 
then obtain g n = Q.92GeV~ 2 and r s = g s j 1 g n = 0.727 from the requirements 

= lAOMeV and Mn = 940MeV, respectively, and the resulting scalar 
diquark mass becomes Md = 600MeV. As we explained earlier, in the TR 
scheme we use the same value = 600Me^, and rewrite the t-matrices r n 
and t s , which are needed to calculate the distribution functions, in terms of 

and Mr, such that they become independent of g n and r s . 

Our results for the valence and sea quark distributions are shown in figs. 
4-7 both for the LB and the TR cut-off scheme. As we have explained 
in sect.l, in order to make contact to the empirical distributions extracted 



10 The current quark mass obtained from the gap equation is m = 5.96MeV in the LB 
scheme. 
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from the measured structure functions, we have to evolve our calculated 
distributions from the low energy scale /i 2 = Q$ to the value fi 2 = Q 2 where 
empirical parametrizations are available. For this Q 2 evolution we use the 
computer code of ref. [[HJ to solve the DGLAP equation in the next-to- 
leading order. (For the Q 2 evolution we use Nf = 3, Aqcd = 250MeV.) We 



will compare our evolved distributions to the parametrizations of ref. |33| for 
Q 2 = AGeV 2 . Both the calculated and the empirical distribution functions 
refer to the MS renormalization and factorization scheme. The value of Qq 
is treated as a free parameter which is determined so as to reproduce the 
overall features of the empirical valence quark distributions at Q 2 = AGeV 2 . 
In this way we obtain a value of Ql = 0.16GeV 2 , i.e; Qo is equal to our 
constituent quark mass Mq. 

Let us first discuss the valence quark distributions shown in figs. 4 and 
5. Although we do not show the results of the pure valence quark model (no 
pions), we note that the input distributions at fi 2 = Qq shown here are softer 
than in the case without pionic cloud effects, that is, the pionic effects reduce 
the peak heights of the valence quark distributions and increase their support 
at low x. The integral over the input distributions shows that at fi 2 = Qq the 
valence quarks carry 92 % (87 %) of the nucleon's LC momentum for the case 
of the LB (TR) cut-off. The rest is carried by the sea quarks. This reduction 
of the peak heights due to pionic effects has a beneficial effect on the overall 
behaviour of the valence quark distributions, although it is insufficient in 
particular for the d quark in the LB scheme. The input distributions are still 
rather stiff even when pionic effects are taken into account, which necessitates 
the use of a low value of Qq in order to approach the empirical distributions 
via the Q 2 evolution. We can expect some improvement concerning this point 
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in a full Faddeev calculation, since in the present static approximation the 
size of the nucleon is too small, corresponding to momentum distributions 
which are too stiff. 

In the LB cut-off scheme, the input distributions are zero for large (and 
also very small) values of x, and therefore the output distributions (at Q 2 = 
AGeV 2 ) show a too strong variation with x compared to the empirical ones. 
On the contrary, for the TR cut-off the input distributions are non-zero in 
the whole region of x, which leads to a smoother behaviour of the output 
distributions, [j This is the same feature as noted in ref. [15| for the quark 
distribution in the pion, and indicates that for phenomenological applications 
the TR cut-off seems to be superior over the LB cut-off. On the other hand, 
as we have explained earlier, the shortcoming of the TR cut-off scheme is that 
the diquark mass must be treated as a free parameter since in this scheme 
mass renormalizations are necessary in order to get finite results. 

In our calculation, the difference between the valence u and d quark 
distributions reflects the scalar diquark correlations in the proton: Since 
the d quark appears inside the diquark and not as a spectator quark (see 
eq. (|3 .£] )), its distribution is given by the convolution of two distributions 
(eq. Q3.9Q ), which is more concentrated at low values of x compared to the 
spectator quark distribution. This is in agreement with the behaviour shown 
by the empirical distributions, and this observation was in fact one of the 
motivations to introduce diquark degrees of freedom also into the bag model 
description of the nucleon structure functions |3J . 



11 For the TR cut-off, the input distributions show a sharp increase when x becomes 
very close to 1. Since the computer code used for the Q 2 evolution |y] requires an input 
distribution which vanishes for x = 1, we artificially modified it for x very close to 1 such 
that it goes like (1 — x) n with some power n. (n — 10 was used in the actual calculation.) 
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We now turn to the antiquark distributions shown in figs. 6 and 7. As 
in the case of the valence quark distributions, the TR cut-off scheme leads 
to an overall better agreement with the empirical distributions than the LB 
scheme. The enhancement of d over u is clearly seen both in the input and 
the output distributions. Since the numerical value of the probability of the 
quark state without pion cloud is Zq = 0.84 (Zq = 0.83) for the LB (TR) 
cut-off scheme, the Gottfried sum ( |3.30|) becomes Sg = 0.262 (Sg = 0.257), 
compared to the experimental value Q reported by the NM collaboration [|TJ 
Sq = 0.235 ± 0.026. Our results for the difference j-^ip ~ fu/p are shown 
in fig. 8. We see that with our value of Ql which has been chosen such 
as to reproduce the overall behaviour of the valence quark distributions, the 
calculated difference is smaller than the empirical one for intermediate values 
of x, but larger for small x. Concerning the ratio u/d, there are also data 



from Drell- Yan processes [35) which give u/d = 0.51±0.04±0.05 at x — 0.18, 
compared to our calculated value of 0.68 (0.70) at Q 2 = 4GeV 2 for the LB 
(TR) cut-off. This, too, shows that the observed flavor asymmetry of the 
Dirac sea is larger than our calculated one in this range of x. 

5 Summary and outlook 

In this paper we used the framework of the relativistic Faddeev equation in 
the NJL model to calculate the quark LC momentum distributions in the 
nucleon. As a first step towards a full Faddeev calculation, we used the 
nucleon vertex functions obtained in the simple static approximation to the 
Faddeev kernel, and included pionic cloud effect approximately using the 



12 It has been shown that the Gottfried sum is almost unchanged by the Q 2 evolution 
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familiar convolution formalism. We can summarize our results as follows: 
First, we have shown a method based on a straight forward Feynman dia- 
gram evaluation, which we believe to be best suited for the calculation of the 
distribution functions in the Faddeev framework. Besides being simple and 
straight forward, the method also guarantees the validity of the number and 
momentum sum rules from the outset. Second, we discussed two regulariza- 
tion schemes which can be formulated in terms of LC coordinates, and which 
preserve the number and momentum sum rules. Based on our numerical 
results we concluded that for the description of the momentum distributions 
in the nucleon the transverse momentum cut-off scheme is superior over the 
invariant mass regularization scheme, which is similar to the situation found 
previously for the momentum distribution in the pion. Third, we have shown 
that the resulting distribution functions reproduce the overall behaviours of 
the empirical ones if the low energy scale for the Q 2 evolution is taken to be 
about the same as the constituent quark mass (400MeV in our calculation). 
Such a low value is required since the input valence quark distributions cal- 
culated in our model are rather stiff, although the pionic cloud effects served 
to soften them compared to the pure valence quark results. In this respect, a 
full Faddeev calculation which gives larger and more realistic nucleon radii, 
as well as the inclusion of higher mass diquark channels (axial vector diquark 
channel) are expected to improve the situation. We have also shown that the 
value for the Gottfried sum obtained in this simple calculation is in basic 
agreement with the experimental one. 

The formulation and results of this work can be used as a basis for at 
least the following three extensions: First, one could easily use the present 
framework to calculate the quark spin distributions, provided that the axial 
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vector diquark ( J = 1) channel is also taken into account. Second, one should 
use the full Faddeev vertex functions to calculate the quark momentum dis- 
tributions. As we have noted in sect.l, the most convenient way for this 
purpose might be first to go to the moment space and then to reconstruct 
the distribution functions. The third extension concerns the case of finite 
density: For a finite density calculation the full Faddeev framework seems 
to be intractable, and approximations like the static approximation used in 
this paper might be unavoidable. For this purpose, however, it is necessary 
first to construct an equation of state for nuclear matter based on the Fad- 
deev (quark-diquark) picture of the single nucleon, similar to the Guichon 
equation of state which is based on the MIT bag picture of the single 
nucleon. The construction of such an equation of state and its applications 



are now under consideration p7]. 
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Figure captions 



1. Graphical representation of the quark LC momentum distribution in 
the Faddeev framework. The single (double) line denotes the quark 
propagator (scalar diquark t-matrix), the hatched circle stands for the 
quark- diquark vertex function, and the operator insertion denoted by 
a cross stands for 7 + 5(/c_ — ± r 2 )/2 for the U(D) quark distri- 
bution. The second diagram stands symbolically for those 2 diagrams 
obtained by inserting the cross into both particle lines in the diquark. 
The quark-diquark vertex contains the isospin operator r 2 . The dia- 
gram where the operator insertion is made on the exchanged quark is 
not shown here since it does not contribute in the static approximation. 

2. Feynman diagrams which have to be evaluated in addition to those 
shown in fig.l due to the presence of the pion cloud. Here the dashed 
line indicates the qq t-matrix in the pionic channel, and the other lines 
are as in fig.l. The second diagram stands symbolically for those 2 
diagrams obtained by inserting the cross into both the q and q lines in 
the pion. The dots indicate all remaining diagrams where the pion line 
is attached to a quark in the diquark. 

3. Graphical representation of the quark distribution within an on-shell 
parent quark. The meaning of the lines is as in figs.l and 2. The quark 
spinor Tq(p) is associated with the incoming and outgoing quark lines 
in all diagrams. The operator insertion in the first diagram stands 
for 7 + <5(p_ — p-x)(l ± t z )/2, and in the other diagrams for 7 + <5(/c_ — 
P-x)(l ± t z )/2 for the u(d) quark distributions. The third diagram 
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stands symbolically for those 2 diagrams obtained by inserting the cross 
into both the q and q lines in the pion. 

4. LC momentum distributions of the valence w-quark in the proton, using 
the LB (solid lines) and TR (dashed lines) regularization scheme. The 
lines associated with Ql = 0.16GeV 2 show the NJL model results, and 
the lines associated with Ql = 0.16GeV 2 show the results obtained by 
the QCD evolution in next-to- leading order from = 0.16GeV 2 to 
Q2 = ^GeV 2 , using A QCD = 0.25GeV and N f = 3. The dotted line 
shows the parametrization for Q 2 = AGeV 2 obtained in ref. |33|] by 
analyzing the experimental data. 

5. LC momentum distributions of the valence ci-quark in the proton. For 
explanation of the lines, see the caption to fig. 4. 

6. LC momentum distributions of the M-quark in the proton. For expla- 
nation of the lines, see the caption to fig. 4. 

7. LC momentum distributions of the d-quark in the proton. For expla- 
nation of the lines, see the caption to fig. 4. 

8. The difference of d and w-quark momentum distributions in the proton. 
For explanation of the lines, see the caption to fig. 4. 
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